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A CLARIFICATION WITH REGARD TO MANUSCRIPT REPORT SERIES NO. 55 
Page 21. (Line 10 from bottom) 
Replace Mo with So 


Page 22. (Line 2 in legend) 
Replace My with So 


Page 23. (The constants in the upper figure should read) 
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Page 80. (Line 15 from bottom and Line 9 from bottom) 
Replace 5 by 6 


Page 89. (Line 10 from bottom) 
Replace 5 by 6. 
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1. AVANT-PROPOS 


Une carte cotidale est une extrapolation géographique de la phase et de 
l'amplitude d'une onde composante, obtenues par ]'analyse harmonique de la 
marée observée en divers points de 1'océan et consiste de courbes de phase 
et d'amplitude constantes. Elle sert a4 interpréter les données disponibles, 
a verifier les résultats d'analyses subséquentes et a comprendre la dynamique 
des marées. Nous pouvons construire des cartes cotidales adéquates pour les 
eaux adjacentes au Canada que 1'on peut interpréter a 1'aide de la théorie 
élémentaire du mouvement ondulatoire dans un océan tournant. La théorie 
suggére 1'existence de structures typiques telles que points d'amphidromie 
et ventres que l'on devrait retrouver dans la nature. Un jeu de cartes co- 
tidales des ondes majeures peut servir aussi a préparer des prédictions 
approximatives de la marée pour tous les points des océans qu'il couvre. 


On a accumulé systématiquement des données sur la marée au Canada depuis 
la fondation du Bureau des Levées de Marées par W.B. Dawson en 1894. 
L'Amirauté Anglaise avait fait auparavant des levées sporadiques dans la 
Baie de Fundy, le Golfe du Saint-Laurent, autour de Terre-Neuve, le long de 
notre cote du Pacifique et dans l'‘arctique. Depuis 1894 on a accumulé 
systématiquement les données et elles ont formé la base d'un grand nombre de 
travaux scientifiques et pratiques. Les cartes cotidales que je présente 
ici constituent une distillation et une interprétation du matériel trés vaste 
qui est maintenant a notre disposition. 


FOREWORD 


Cotidal charts constitute a geographical extrapolation of the major har- 
monic constants obtained from a harmonic analysis of the tide observed at 
points in and around the ocean; they consist of contours of equal phase and 
equal amplitude. They are of value in interpreting the data available, for 
checking the results of further analyses and for understanding the dynamics 
of tidal motion. It is possible to draw fair cotidal charts for the waters 
adjacent to Canada which may be interpreted using the results of the theory 
of elementary wave motion in a rotating ocean. Such theory suggests struc- 
tures which should be found in nature such as amphidromies and antinodes. A 
set of cotidal charts of the major harmonic constituents can also be used to 


make approximate predictions of the tide at any point of the oceans which 
they cover. 


Tidal information about Canadian waters has been accumulated systematical- 
ly Since the establishment of the Tidal Survey with W.B. Dawson in 1894. 
Previously the British Admiralty had made sporadic measurements in the Bay 
of Fundy, the Gulf of St. Lawrence, around Newfoundland, on the Pacific west 
coast and in the arctic. Since 1894 the information has been systematically 
compiled and scrutinized to become the source material for numerous scientif- 
1c and practical investigations. The cotidal charts which I present here 


constitute One such result; a distillation and interpretation of the very 
extensive material now available. 


2. COTIDAL CHARTS 


Cotidal charts provide a visual summary of available information about 
the major harmonic constituents of the vertical tide. If a plot is made on 
a geographic chart of the amplitude and phase of these constituents at the 
stations available, it is readily apparent that they vary smoothly from point 
to point and it is possible by simple extrapolation to draw contours over 
the water body. No theory is involved at this stage. 


Elementary physical theory indicates that the tide is a long wave phen- 
omenon. Harmonic constituents have a precise frequency and a basin of a 
given configuration will respond in a definite way to a forced oscillation: 
nodes and antinodes will be formed at specific locations. These should be- 
come apparent during the inspection of the data. The existence of a node is 
revealed by a rapid change of phase around a given point and a marked de- 
crease in the amplitude there. An antinode is found in an area where the 
phase is virtually constant while the amplitude passes through a maximum. 
The end product of this task of interpretation, interpolation and extrapol- 
ation is a cotidal chart. The data are seldom dense enough to draw it unam- 
biguously so that a cotidal chart contains a certain element of arbitrariness. 
Nevertheless, cotidal charts are useful tools in study of the tide and one 
should never hesitate the prepare one even if the information available 
appears sadly lacking. 


In briert: 


1) Cotidal charts help summarize the information presently available 
on the tide in a given area. 


2) They may be used to check the results of analyses of short sets 
of observations at temporary gauging sites; gross errors in the recordings 
or in time keeping may be picked up by plotting the analyzed values on the 
appropriate charts. They also form the basis for more precise and plausible 
sets of cotidal charts when additional information becomes available. 


3) They help in pointing areas where the tidal coverage is inade- 
quate and where further exploration is necessary. Canadian waters are 
extensive and only by careful planning and systematic execution will it ever 
be possible to obtain a general understanding of the tide in our waters. 


4) Wherever the tide has a good semidiurnal character (this will 
be defined more precisely further on), the Mz chart alone suffices to give a 
good estimate of the mean amplitude of the tide at a given site (the range 
is twice the amplitude) and of the mean time of high and low water. 


5) A correctly drawn cotidal chart reveals the structure of the 
tidal wave of the given frequency in the basin under consideration and helps 
stimulate a search for a theoretical confirmation of the assumed structures. 


6) Cotidal charts may be used to predict the tide in the particular 
area. 


7) Cotidal charts help control the results of numerical modelling 
of the tide in the ocean. 


2.1 Geographic chart 


I have used a Miller projection as a background geographic chart 
(frontispiece). This is an intermediate cylindrical projection. The merid- 
ians as well as the parallels are straight parallel lines which have been 
devised for better representation of high latitude regions. The parallels 
are spaced much closer than on Mercator's to reduce excessive latitudinal 
expansion, and farther apart than on the Lambert Equal-Area Cylindrical, to 
avoid extreme polar flattening. 


No projection can display perfectly on a plane a large portion of the 
earth's sphere; the Miller projection in particular gives more emphasis to 
northern areas, and since we have much more exploratory work to do in our 
northern regions its use helps emphasize where more measurements are required. 
In any case the southern regions where more measurements are available are 
presented here in darger scale. 


3. HARMONIC CONSTITUENTS: AMPLITUDE AND PHASE 
The vertical tide z(t) may be adequately represented by the expression: 


n 
ie} eS + YF : -UC-d. 
Zt) Ze a A. Cos (ot a;) (1) 


where Z,) is a constant which denotes the reference level, which in ocean- 
ography represents the mean sea level (whenever this quantity can be adequate- 
ly defined and calculated); on navigational charts Zg is an arbitrary 
elevation chosen in such a way that the observed water level seldom falls 
below it. The second term is a superposition of n pure harmonics of ampli- 
tude A, frequency o and phase a; it represents the tide proper. (1) In- 
dicates that at most places the tide may be represented by a sum of sinusoids 
of known frequencies, but whose amplitudes and phases are characteristic of 
the site. A harmonic analysis determines the amplitude and the phase of 
each harmonic component (constituent) from observations on the local varia- 
tions of water level. The number n of such constituents is very large (over 
500) ; in practice a handful of harmonics in the diurnal and semidiurnal 

bands is much larger than the others combined and suffice to give a good 
approximation to the tide. These are given in Table 1. Other constituents 
such aS U2, V and L2 also contribute significantly to the tide, but as their 
resolution and interpretation are more difficult we exclude them from the 

set of charts we wish to present. We also exclude the charts of P] and Ko, 
because these constituents have frequencies close to Kj and So respectively. 


ee contribution of a given constituent to the local tide is, as already 
said, a pure sinusoid of amplitude A and phase a which we illustrate in 
Figure 1. We have found it preferable to replace the time scale t by the 
phase change ot which enters in (1). We notice immediately that the phase a 
will vary with each time origin th chosen for a specific analysis. This 1s 
a source of confusion and an accepted convention in tidal analysis is to 


Table 1. Frequency of the major tidal constituents. 


Diurnal band Semidiurnal band 
Constituent Frequency Constituent Frequency 
(deg/h) (deg/h) 
0, 13°94 No 28.44 
Py 14.96 My 28.98 
K, 15.04 So 30.00 
K., 30.08 


Fig. 1. Relation between the local constituent A cos (ot-a) and the 
corresponding component of the tidal force A' cos V(t). 


replace the phase lag a by the Greenwich phase lag g which is defined in the 
following way. The phase of the component of the tidal force corresponding 
to the specific constituent can be calculated from known astronomical data; 
it is called V(t), the “astronomical argument" and it has value V(to) at the 
chosen time origin for the analysis. V(to) is calculated routinely in the 
course of a harmonic analysis. The Greenwich phase lag g is defined by the 
expression: 


Sr aaa oes (2) 
or better; 


g =V(t.) +a (3) 


Figure 1 illustrates the meaning of g: it is the phase lag between the in- 
stant the component in the tidal force reaches its peak and the instant of 
maximum contribution of the constituent to the local vertical tide. Physical- 
ly it is the delay between the time a force is applied and the time it causes 
maximum deformation. The phase lag g has the same value for any choice of 
the time origin ty. However the value of g also depends on the time zone Z 
used to determine the time of the culmination of the component tidal dis- 
placement; in order to resolve this ambiguity, all the phase lags g in the 
charts presented have been calculated for the Greenwich time zone 


Z=0, 


so that the g's contoured in the charts are directly comparable from one site 
to the other. The translation of the Greenwich phase lag for a local time 
zone Z=r is given by the formula: 


da Gerla (4) 
As an example the Greenwich phase lag for Ky al ViGtoni ads: 
g = 269.2° 


Since Victoria lies in time zone Z = +8, the Greenwich phase lag for Ky 
(frequency 15.04°/h) in local time is: 


Jo = 269.2 - 8x15.04 = 148.99 


The amplitude A of the constituent is independent of the time scale and of 
70 the level of reference used. The value A and g of the constituents at 
€ gauging stations where the tide was observed form the source material of 


the cotidal charts. 


4, PERIODIC WAVE MOTION 


It is useful in the preparation of cotidal charts to have some under- 
Standing of periodic wave motion in a basin, which understanding, in turn, 
can only be obtained with the help of a formal solution of the linearized 
equations of hydrodynamics. In rectangular components these have the form, 
taking account of the rotation of the earth and of friction: 


ou OZ 
ey = =U (5) 
ot OX 
OV oZ 
etd PY (6) 
ot oy 
du OV oZ 
H{—+—]+—=0 (7) 
OX oy ot 


where u, v are the east and north components of the current vectors; Z is 

the vertical displacement of the surface; f = 22 sin (in which 2 is the 
speed of rotation of the earth and ¢ the latitude) is the Coriolis coeffic- 
jent; g is the acceleration due to gravity; H the depth (assumed constant) 
and r is the linearized coefficient of friction. Motion is presumed two- 
dimensional and homogeneous in the vertical (barotropic); (x, y) denotes the 
coordinates of the point of measurement in flat space and t the time variable. 
For periodic motion (the case of the tide), the time variation of u, v and Z 
can be represented by: 


10ot 
Uae Via Oke 


so that 2+ becomes io and (5), (6) and (7) take the form: 


ot 
OZ 
(io +r) u - fv = -g — (8) 
OX 
OZ 
fu + (io + r) v = -g — (9) 
oy 


ou OV 
hg Se ss fo (10) 
OX ay 
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(8) to (10) form a set of linear first-order differential equations in the 
unknowns u, v and Z defined over a basin whose dimensions are specified for 
each problem to be solved. In the present study we presume that the basin 

is rectangular. Their solution in a finite domain consists of Kelvin and 
Poincaré waves that can be combined linearly in order to satisfy specific 
boundary conditions. In spite of the ready availability of complete and 
exact solutions of (8), (9) and (10), they are so involved that it is prefer- 
able to work toward them in stages of increasing physical complication. I 
begin with one-dimensional motion. 


4.1 One dimensional motion (no rotation) 
4.16). No Frvecior 


Equations (8), (9) and (10) reduce to: 


OZ 
jou = -g — ‘aka 
OX 
ou 
H =—= sioZ (12) 
OX 


where we made v = f = r = O and the basin becomes a canal. The solution of 
(11), (12) with the time variation is: 


7 = (Ae 1K me B e” 1KX) lot (13) 


g : ; , 
4° (Ae 1KX _B e7 1KXy aict (14) 


where cs the wave number, stands for o/YgH. A and B are arbitrary constants 
of integration to be determined from the boundary conditions. The physically 
Significant part of the solution is the real part of the expressions for Z 
and u. We note that: 


= 
II 


e “=8c0s a 4 7 san a (15) 


where a is usually expressed in radians and is dimensionless. 


4.1.1.1 Progressive waves 
If the boundary condition is such that: 


A=) or B= 0 


the solution (13), (14) reduces to: 


_= B oi (ot-kx) x 5 


el 
bd os oye oi (ot-kx) ings ay’ oi (ottkx) 
Ih H 
H 


which, in real form, assuming A and B real, become: 


_ A oi (ottkx) 


Z, = B cos (ot-kx) (16) Z, = A cos (ottkx) (18) 
g g 

u, = B F cos (ot-kx) (179 Up = -AQ— cos (ot+kx) (19) 
H 


These solutions represent oscillatory waves travelling in the direction of 
increasing or decreasing x. In order to visualize this we study the profile 
of the elevation and the current along the canal (Figure 2a) between kx = 0 
and kx = 360° for various time steps for the solution (16) and (17). The 
heavy line in the surface represents its displacement with respect to its 
equilibrium position. The arrows denote the direction and intensity of the 
associated current; a dot indicates a position of null current. The four 
consecutive figures for times corresponding to 0° to 270° of phase indicate 
that the surface appears to be affected by a wave travelling from left to 
right. The associated currents are null at the nodes, positive and maxima 
at the crests, and negative and minima at the troughs. The second solution 
would represent a wave travelling from right to left with the same velocity. 


The wavelength L is obtained from: 


Z(x, +L) = Z(x)) or cos (at-k(x,#1) =1 COS (ot-kx,) 


which implies 
kKLe=*29 or L=— 
The period T of the wave is determined from the condition: 


Z(t+T) = Z(t) or cos (o[t+T]-kx) = cos (ot-kx) 
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5 t=0° 


0° 90° 180° ZO? 360° 


1) =a cos (Kx -gt) 


g 
u = [=] a cos (Kx —gt) 


o° 90° 180° 270° 360 K= 


ot =|80° 


ot=270° 


‘OF 90° 180° 2102 360° 


Fig. 2a. Profile of the surface of a fluid in a canal and corresponding 
horizontal velocities at various times when it is disturbed by a travelling 
wave. 
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which implies 


Q 


Finally the velocity c of displacement of the crest is obtained from the 
condition: 


— (ot-kx) = 0 or o=k—=kc 


Therefore: 


The second solution gives a wave with the same wavelength, period and 
velocity, but with displacements in the opposite direction. 


4.1.1.2 Standing waves 


PEehc= bore? A= —b. theesolution (13) (14) takes ‘the 
special form: 


lL = hCos bemos ot (20) oe Asam «x COS. OL (22) 


g 
Kh yen kx cos (ot-1/2) (21) U» 
H 


where BO = 2K Ave =e Dad 


A" cos kx cos (ot+1/2) (23) 


Los 
al 
II 


These two special solutions represent standing waves. In (20) (21) the 
current u vanishes at: 


x=10) 


This point can be considered as the position of a wall against which a wave 
travelling from the right is perfectly reflected to return to the right with- 
out any change in amplitude. This interpretation can be supported algebrai- 
cally noting that: 


14 


iy = A' cos kx cos ot = A cos (ot+kx) + A cos (ot-kx) = 2A cos ot cos kx 
ee — ee 
wave travelling wave travelling 
from the right to the right 


The second solution also represents a standing wave with perfect reflexion 
at kx = 71/2. In a standing wave the water level oscillates vertically with- 
out apparent displacement of the wave crest. The associated current is 
maximum at the points where the amplitude of the vertical displacement is 
zero and vice versa. The wavelength is still: 


L =. 20/k and the period is still T= 21/0. 


but the velocity of displacement is 0 (Figure 2b). Standing and travelling 
waves represent extreme cases of periodic motion in one dimension. Any 
intermediate type of motion can be considered as a superposition of two 
travelling or two standing waves with different amplitude and phase. 


Tidal motion in the ocean and in neighbouring basins has a standing wave 
character and we will be able to interpret the cotidal charts by representing 
the motion by a solution of the form (20) (21) modified for the effect of the 
earth's rotation and of friction. 


4.1.2 Wavelengths of tidal waves 


The cotidal charts are presented for the constituents 01, 
eae i and So. Their frequencies and associated periods are as shown 
able 2). 


Table 2. Periods of the constituents considered 
in the cotidal charts. 


Constituent Frequency Period 
fe} Hi 
wh radians/sec h 
0, 13.9430 9 © 67e08" 25.8 
K, 15.0411 .729x1074 23.9 
No 28.4397 1.379x1074 12.7 
M, 28.9841 1.405x1074 12.4 


- 30.0000 1.454x1074 12.0 


Lbs 


Kx 


o° 90° 180° 270° 360° 


y - acos Kx cosot 


et a a sin Kx sinot 


oO 
Keo 
{gH 
ot=|80° 
Ed & 
T= 
oO 
° ° ° fe} 
(6) 90 180 270 360° : oT 
nik 
a : 1@) 
6ot=270°3 et 
_ a. Kx 
Ops 90° 180° 2002 360° 


Fig. 2b. Profile of the surface of a fluid in a canal and corresponding 
horizontal velocities at various times when it is disturbed by a standing 
wave. 
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We calculate the wavelengths L and velocities c of the constituents M2 and 
Kk] which predominate in the semidiurnal and diurnal bands for various values 
of the depth H (Table 3). Although tidal motion behaves as a standing wave, 
in general we still calculate the velocity c of a progressive wave. The 
reason is that the tide may assume the character of a travelling wave in 
narrow and shallow rivers; the effect of friction is such that the velocity 
of displacement of the tidal wave is less than the value indicated in the 


table. 


The order of magnitude of the wavelength is more relevant to the inter- 
pretation of cotidal charts. The wavelength of the diurnal tide is twice 
that of the semidiurnal, so that the distance between diurnal nodes should 
be twice as large as between semidiurnal nodes. If we find a node, semi- 
diurnal or diurnal, it provides a clue as to the location of the node of the 
other species. The distance between nodes of the same species should be 
half a wavelength (for a constant depth). 


Table 3. Velocity and wavelength of the tidal waves 
associated with the constituents My and Ky. 


Depth Associated velocity Mo Ky 
for<a 
free travelling wave Wavelength Wavelength 
H c= Noi L=2n/k=Tc L=Tc 
m km/h km km 
10 36 447 862 
50 80 994 1915 
100 TT3 1404 £705 
500 ee 3130 6031 
1000 Sor 4435 8545 
5000 797 9901 19076 


4.1.3 Construction of cotidal charts from mathematical models 


We wish to utilize the elementary theory developed ab 

draw cotidal charts corresponding to travelling ee Birerte ene Heelan 
tidal chart consists of the field of amplitudes and phase lags in the time 
dependent part of the solution, whether it be the elevation Z or the current 
w In the solution for a travelling wave, i.e. (16) and (17), both Z and u 
ae the same time dependence, kx being the phase lag; their amplitude differs 
a a ere factor Vg/H. Here one cotidal chart suffices for both Z andu. 

€ solution for a standing wave, (20) and (21), gives no phase lag for Z 


w@ 


but the amplitude is affected by the factor cos kx; this factor determines 
areas where the amplitude is positive or negative. Since we keep the ampli- 
tude always positive in our charts, we treat a negative amplitude as con- 
sisting of a positive amplitude with a phase lag of 180°, j.e.: 


-A cos ot = A‘cos (ct-180°) 


In a canal, the point where cos kx = 0 determines the node which takes the 
form of a line across it. The solution for the current indicates a constant 
phase lag of 90° with respect to the vertical tide, but the possible sign 
change of sin kx also brings an additional phase lag of 180° in the zones 
where sin kx is negative. 


With these remarks we go on to search for a solution of a standing wave 
with friction and then will draw cotidal charts for all the cases studied in 
order to see their appearance and their modification by the introduction of 
additional effects. 


4.1.4 One-dimensional standing wave motion with friction 


In the case of one-dimensional motion with friction 
(v = f = 0) the equations have the form: 


OZ 
(iotr) u = -g — (25) 
OX 
ou 
Hts SE iG7 (26) 
OX 


and the solution has the same form as (13) and (14) except that k is replaced 
by k' where: 
Oo 
Kk. == yl -ir/o 
VgH 


The form of this factor indicates that it is no longer an elementary matter 
to separate the real and imaginary terms in the solution; it is useful to 
retain the time dependence of the time in the form ei°t until the very last 
Stages of the calculations. Whatever complex factor is present in the sol- 
ution it can always be written in form Re~'", where r is an added contrib- 
ution to the phase lag. 


The solution of (25) (26), which represents a standing wave with perfect 
reflexion at the origin (u=0 at x =0), is: 
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A cosjuk.x.COS, ot (27) 


N 
Il 


g sin k'x cos (ot-1/2r) 
u = Aa, ——._ 
H yl-irso 


where we have dropped the prime in A'. Since k'x is a complex number, its 
cosine or sine consists in fact of products of trigonometric and hyperbolic 
functions. It follows then that the constant A is complex in general and 
this must be taken into account in the calculation of the effective phase 
lags of the elevation or of the current. The current satisfies the condition 
of perfect reflexion at x = 0 because sin k'x = 0 for x = 0 even if k' is 
complex. On the other hand cos k'x has no zero whatsoever, which in turn 
indicates that there is no node possible for such type of motion. We call 
it a standing wave with perfect reflexion; it could be considered as well 
the superposition of two waves travelling in opposite directions, with the 
wave reflected from the wall having a smaller amplitude. We now draw co- 
tidal charts for a travelling wave and a standing wave without friction and 
then another set for a standing wave with friction. We show only the chart 
for the elevation of a travelling wave. 


The cotidal chart for the travelling wave (Figure 3a) shows a canal open 
at the origin; it contains only cophase lines which increase for increasing 
x. There are no coamplitude lines because the amplitude of the wave is 
uniform everywhere. In the second pair (Fig. 3b) of cotidal charts the 
canal is closed at x = 0. The charts are those of the vertical tide Z and 
the associated horizontal current u in a standing wave. Cophase and co- 
amplitude lines coincide in many points; in this case the coamplitude is 
pulled out of the canal and the corresponding amplitude is written outside 
as well. In the vertical tide we get a maximum amplitude at the origin and 
at the position kx = 180°. The node is at kx = 90° and it is marked by a 
coamplitude line of 0. The corresponding cotidal chart of the current shows 
that it has uniform phase between x = 0 and kx = 180° where the node of the 
current is to be found. It has maximum amplitude at the position of the node 
of the vertical tide. The other pair (Fig. 3c) of charts show the same 
Standing wave, but this time affected by friction. The chart for the verti- 
cal tide shows that the node at kx = 90° has disappeared and has been re- 
placed by a minimum of .2. The jump of 180° at the node is now replaced by 
a rapid change of 340 to 220°. Similarly for the current, the abrupt change 
in phase from 90° to 270° is replaced by a change from 70 to 10° in the vic- 
inity of the node of the current. The latter charts are more physically 
reasonable. They show half a wavelength in a standing oscillation. Numerous 
observations are available on Z which permit us to verify such an hypothesis; 
whereas the few data on u forbid the elaboration of similar cotidal charts 
for the horizontal tide. 


4.2 Two-dimensional motion with rotation and friction 


The equations to solve are (8) to (10) and their general solution has 
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(a) (b) (c) 
Z 7 p s ’ 
—Coamplitude 
= Spe ae 
Ky 
— OPE OSED 
N cL 
Travelling Wave standing ee 4 ee ee 


NO FRICTION FRICTION 


Fig. 3. Cotidal charts in a canal half a wavelength long (one-dimensional 
motion). Continuous line - cophase contour. Dashed line - coamplitude 
contour. 


a) Travelling wave: vertical displacement Z. 

b) Standing wave, no friction: vertical displacement Z and 
horizontal velocity u. 

c) Standing wave affected by friction: Z and u. 
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the form: 
Yi (Ae i KIX 4: Bew KIX) (ca! Key + De 1K2Y) a (29) 
where 
» (o-ir)*-F 
kK, a ks = 
gH 


ky, Or ko, usually Torms ad set of eigenvalues whenever a boundary condition 
of zero velocity is imposed on the x or y boundary; the corresponding set of 
solutions is known as Poincaré waves. If ky (or ko) is pure imaginary (this 
is determined by the width and depth of the sea), the Poincaré waves are 
damped rapidly away from the boundaries and are not significant over the 
body of the sea. 


In addition to (29), the equations (8) to (10) also have a solution which 
satisfies: 


< 
Wt 
S 


throughout the domain; this is the Kelvin solution. It is the equivalent of 
that which we have been studying in one-dimension, but with the additional 
effect of rotation. It represents a very simple mode of motion and as we 
shall see, seems to be preferred by the ocean. 


The Kelvin solution for (8) to (10) analogous to the one-dimensional 
solution (27) (28). qs: 


jot 


N 
II 


A cos (k'x-im'y)e 


g ] : 
=a sin (k'x-im' i(ot-1/2) 
H Yl-ir/o co ae 3) 


O f ] 
where k! = yl -ir/o i 4 
: = —_—————-_ and where 
V9H {GH V1-ir/o 


m' represents the distorting effect of rotation. The prime in both parameters 
k and m indicates that linearized friction is also taken into account. 


: We give (Fig. 4) schematic cotidal charts of the solutions (30), and (31) 
or the particular case of no friction (r = 0); they should be studied in 
comparison with those of one-dimensional standing wave with no friction. 
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Fig. 4. Schematic cotidal charts for a standing Kelvin wave (no 
friction); they are the rotary counterpart of chart (b) in Fig. 3. 
The nodal line has been transformed into a point (amphidromy) by 
the rotation and the cophase lines now radiate from it. The node 
of the current (at X1) corresponds to the antinode of the vertical 
displacement and vice versa. 


The new chart of Z may be compared with the original one being rotated 
counterclockwise and we see that the nodal line is now a point (amphidromy). 
Along the central line of the basin the lines follow exactly the same 
progression as in one dimension. The same comparison can be made with the 
charts for the current u. We do not draw the charts as far as the boundary 
x = 0 in two-dimensions because the Kelvin solution is not the only one 
present near the boundary; an infinity of Poincaré waves aiso exists in that 
vicinity in order to ensure the satisfaction of the condition of zero flow 
all across the boundary at x = 0. 


The diagrams given are valid for the northern hemisphere; the crest of 
high water rotates counterclockwise around the node. We notice that maximum 
currents are still found at the node of the vertical tide and vice versa. 

As already said, we have little opportunity in practice to verify the cotidal 
chart of the current u. 


We give (Fig. 5) the cotidal charts for a standing Kelvin wave for vari- 
ous degrees of friction in a specific basin for the constituent Mo and ill- 
ustrate exclusively the vertical displacement Z given by the solution (30). 
The contours are in arbitrary units and are not directly comparable from the 
case of no friction to those with friction. They indicate qualitatively 
that along the canal the change in amplitude on the "left side" of the canal 
(with respect to the reflecting wall) becomes less marked for increasing 
values of the friction while on the "right side" the amplitudes pass through 
an increasingly deeper and more abrupt minimum. The gradual displacement of 
the node by increasing friction can be confirmed quantitatively by searching 
the position of the node of solution (30). 
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Fig. 5. Cotidal chart of the vertical displacement Z for a standing Kelvin 


Wave corresponding to the frequency of Mo in a basin at latitude 45°N, width 
200 km, depth 100 m, closed at x = 0 and of indefinite length. (a) Without 
friction. (b) With weak linearized friction, r = 10-°/sec. Ccj° (Wich 
moderate friction, r= 5 x 10-9/sec, (d) With strong friction, r = 10-4/sec. 


320° 300° 280° 260° 240° 220° 
160) Km S ——— T ‘ s es = 
pel 
/ fe = bie 
7 4 ees 
“x 3402 : => 
6/ a * 
/ / 5/ 4 Ui — 
/ ae Py 
_— / ! iS Noe 
(c) i ff i . ot Y 
/ | | ae ie 
7 | tf winds if / 
| : 
eae \ 
0° | | | | \ / \ 
| ; \ - 
=f i ake ome | 2 Ne ~| 
me 100 200 #3001 
20° & 100’ 140° 160° 180° 
x — 
NODO A o = 345 Km 


300° 280" 260° 240° 


08 340° 320° 
100 a ca T 
ve < a ek 
/\ 
O F (d) 
Y | 
{ > <a 
-l00 + ~-L- 
100 200s ths ah ye 400 
0° 340° 320° 20° BO. j60° 220° 
SS 
-4 -| 
R= 10 seg NODO A Xo = 356 Km 
Yo = -190 Km 
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The node of the vertical tide Z will be found at the point (X9> Yo) whose 
position is such that cos (k'X9-im'yo) = 0 or such that: 


k'x,-im'y, = (2n+1)-1/2 where W=Osel, 25 Soucy 
For n= 0, the position of the node is given by: 


k'x,-im'y = 42 (32) 


and since k' and m' are complex numbers, (32) is equivalent to two linear 
equations which, when satisfied, give the coordinates of the node. In terms 
of real numbers: 


k' = a-ib in = eid where a, b, c and d are real numbers using 
their definition in (31). Equation (32) is then equivalent to the linear 
equations: 


AX, + oy =e bx, teeyee= 0 


equating real and imaginary parts. Their solution gives the coordinates of 
the node displaced by friction: 


TT C TT b 
SS Ve > 
OY Peach) : 2 (ac-bd) 


Usual values of friction are such that a > b and c > d so that the numerator 
is positive.__Keeping terms to first order with the approximation 

a = k = o/ VgH we get Xq = xg (I+bd/ac), so that the node is pushed away 
from the reflecting surface by friction. We have written xg for the original 
coordinate of the node for the case of no friction: 


yar 


1 
Pome 


X 


\ 
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Similarly the negative sign of Yo indicates that friction pushes the node to 
the right of the reflecting wall (northern hemisphere): 


wall 


ra 


When the value of yo is larger than the half width of the basin, the node 
falls outside the basin and is said to be "degenerate" or "virtual". This 
iS a common occurrence for either strong friction or very narrow channel. 


We will see that with these elementary notions of periodic wave motion, 
of nodes and of standing Kelvin waves modified by friction, we are able to 
interpret intelligently the cotidal charts we have drawn for Canada. 


5. TIDAL PATTERNS 


The tide observed at various harbours has traditionally been described 
as semidiurnal, mixed or diurnal depending on whether its succession of high 
and low waters occurs at regular intervals of six hours, at irregular inter- 
vals (which can stretch between 3 and 12 hours) or at regular intervals of 
12 hours. The type of tide found at a given site depends on the relative 
amplitude and phase of the harmonics making it up. The semidiurnal constit- 
uents overwhelm the diurnal ones in a semidiurnal tide, etc. 


If one scrutinizes the tide at many sites it becomes obvious that there 
is a smooth and gradual transition from one record to the next, so that one 
is faced with the need to develop a criterion which permits the characteriz- 
ation of the tide at a given site and the development of a system of class- 
ification of the various patterns. A criterion much used in the past is the 
quotient of the sum of the amplitudes of Kj and 0] to that of Mo and So, as 
it served an immediate measure of the relative importance of the diurnal and 
semidiurnal signals. The use of such a criterion requires a preliminary 
harmonic analysis. Advances in computer techniques now permit the objective 
description of whatever phenomenon that can be described by numbers with the 
help of so called "patterns" (Batchelor 1972). These can be calculated from 
records and from these values we may draw charts of patterns similar to co- 
tidal charts and describe visually and quantitatively the various types of 
tides encountered (Godin 1977). In our restricted use tidal patterns are 
sets of numbers and labels which help describe, recognize and classify tidal 
records. They are calculated from second moments of first and second differ- 
ences of observed tidal elevations spaced 3 or 12 hours apart. These moments 
take characteristic values for various types of tides and one may put some 
order into them by referring them to predefined "pure" tides. Once the mom- 
ents are calculated and compared with those of the pure tides the pattern of 
any record will fall at some point inside the scale given below: 


io TAS pel N ule a 


where S is semidiurnal, N is normal, M is mixed and D is diurnal. The S and 

M categories are always accompanied by a number smaller than one; the reverse 
holds for N and D. The pattern numbers move smoothly from S to N and M to D; 
we call the points of junction, 1.0N and 1.0D. The boundary between M and N 

is not as sharp; patterns with << 1M could very well fall into >> IN just 
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for a very slight change in the value of: the second moment of the 12 hour 
first difference. This however causes little practical aitetCurty.. |G. 
dicate this situation we mark by a wavy line the approximate zone of demar- 
cation through which the tide goes from a mixed to a normal character and an 
arrow indicates the direction where the normal region lies. 


Note that the classifications S$, N, M and D do not exactly coincide with 
our intuitive conception of these words. First we recall that tidal patterns 
merge smoothly one into the other and that there are no sudden jumps in 
character as the different labels in these categories might suggest. Essen- 
tially one could say that the S, N and M patterns describe "semidiurnal" 
tides as we conceive them intuitively, except that in the S category the 
diurnal inequality is little marked and that in the M category it becomes 
increasingly noticeable. These classifications were elaborated from the 
spectrum of the tidal potential. The N classification contains tides which 
have diurnal and semidiurnal signals in nearly the same proportion as in the 
tidal potential; if less ( <1), it is called S (semidiurnal) because this 
part of the signal is larger than in the applied force. If the diurnal con- 
tribution is more than that in the potential, we enter eventually the M and 
D categories. 


Figure 6 shows the predicted tide for April 1978 for a set of stations 
whose patterns span the scale just indicated. We have chosen stations for 
which we have good records; aS a consequence we could not illustrate all the 
discrete values of the patterns. An inspection of the figure allows us to 
follow the evolution from a highly semidiurnal tide to a very mixed tide with 
prolonged episodes of purely diurnal activity. It must be kept in mind how- 
ever that the profiles we obtain for April will not be exactly reproduced 
during other months. Inspection of the profiles indicates that the success- 
jon of high and low water becomes appreciably distorted beyond patterns 
> .6M. This marks the practical limit beyond which the computation of the 
conventional lunitidal intervals becomes questionable. As well this puts us 
in a position to define precisely what we mean by a "semidiurnal" tide: these 
are the tides which fall in the categories N, S or whose pattern number is 
< .6M. For all the areas of the chart where this criterion applies we need 
to use only the Mo chart to deduce the mean tidal amplitude (and range) as 
well as the mean time interval between the lunar transit and the occurrence 
of high or low water. We use Tofino as an example. Its station number is 
8615 and its geographical location is 49°09'N 125°55'W; the pattern number 
is .161M SO that the tide at Tofino is "semidiurnal", just as defined. The 
Mo chart gives A= 99 cm g = 241°. As a consequence the mean amplitude of 
the tide at Tofino should be 99 cm, the mean range 198 cm, the mean lunitidal 
interval for high and low water should be 8.3 and 14.5 h respectively (the 
low water for Mo follows the high water by 6.2 solar hours). Statistical 
calculations on the observed times and heights of high and low waters at 
Tofino indicate that the mean range is 213 cm while the mean lunitidal in- 
tervals for high and low water are 8.3 and 14.5 h respectively; the 
M2 chart supplies a fair approximation to the mean range and the mean luni- 
tidal intervals wherever the tide is "semidiurnal'. 


ik oe paneer numbers used to elaborate the chart have been derived from 
thea Of analyzed constituents for the given stations; they could have been 
‘ ed as wel] from the hourly observations on the water level. We have 
chosen the first alternative because it minimizes the effect of the noise. 
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Fig. 6. Set of predicted levels for April 1978 at various stations to illus- 
trate some actual tides corresponding to the indicated patterns. The 
stations chosen cover patterns from .4SD (almost pure semidiurnal) to 2.0D 
(almost pure diurnal). (a) 391 SD. Woods Harbour. Good semidiurnal tide 
with a slight diurnal inequality. (b) .483 SD. Halifax. Quite similar to 
Woods Harbour with a more marked diurnal inequality and more intense change 
from springs to neaps. (c) .608 SD. Riviére du Loup. Still a good semi- 
diurnal tide with an even more marked change from springs to neaps. 

(d) 1.004 N. Holy Rood. We enter the normal category where the spectrum 
of the constituents is very similar to that of the tidal potential. (e) .195M. 
Riviére au Renard. We enter the M (mixed) category where the diurnal portion 
of the spectrum of the constituents becomes larger than that in the tidal 
potential. The diurnal inequality is quite marked during the two semimensual 
maxima in the lunar declination which occur around hour 288 and hour 624 of 
the month. (f) .306 M. Pictou. (g) .501 M. Seattle, Washington, USA. 
(h) 602 M.. Lund. -(i) 3s07 M.St@abriola Pass... (7). 1.003'D. Cowichan: 
We enter the D (diurnal) category where the probability increases that there 
will occur a single tide over a span of 24 hours. (k) 1.205 D. Saanichton 
Bay. (2) 1.401 D. Esquimalt. (m) 1.603 D. Caissie Point. This station 
lies to the southwest of the semidiurnal point of amphidromy which exists in 
Northumberland Strait. (n) 2.021 D. Gorge Victoria. Apparently the record 
with the most diurnal character available in our files. Both the diurnal and 
semidiurnal bands of the mixed tide entering from Victoria are damped while 
it propagates into the Gorge but the semidiurnal band is damped more rapidly. 
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Shallow water tides have patterns which fall in the range just described, 
but the fact that their low water is clipped at about the same height for 
all tides necessitates the calculation of an additional index which is de- 


fined by: 
mo(d,)/.252 = | 


where mo(d3) is the second moment of the second differences between hourly 
observations spaced 3 hours apart. The value of I for any shallow water 
stations will be > 1 while for other stations it is < 1 (and not calculated). 


6. COMMENTS ON THE COTIDAL CHARTS 


We have drawn the charts for the constituents M2, S9, No, Ky and 0] in 
order of frequency, as well as charts of the tidal patterns. Techutcaliie we 
should have also drawn the charts of P] and Kg which are also constituents 

of importance. The reason for not doing so is that these two constituents 
are not separable from Kj and So from short records, so that they must be 
inferred. Once their values are laid out on a chart, a good proportion of 
them cannot be treated as true data but only as a direct inference from the 
constituents mentioned. Since on the other hand their value must be includ- 
ed in the preparation of predictions, we may use the ratios and phase dif- 
ferences given (Fig. 7) to calculate them from the amplitude and phase of 

Kj and So. The figure gives the vectors r and ¢, where r is the amplitude 
ratio and $ the phase difference, of P] to Kj or Ko to S2 actually found in 
the course of one year analyses at stations where there was sufficient data 
to separate them. The vector average of all the r and for a given region, 
as seen from the diagram, is sufficient to give an adequate estimate of these 
additional constituents. 


The constituents actually illustrated fall into two species, diurnal and 
semidiurnal. M2 predominates in the semidiurnal band and kj in the diurnal, 
while the maps of Na, So and 07 tend to be quite similar to those of M2 and 
Kj. As a consequence our comments will be limited in general to the maps of 
M2 and K}, although the maps of 0], Na and S? are also necessary to make 
predictions at a given point of the ocean. 


In areas of the maps where the patterns fall to the left of .6M, one may 
use the map of M2 exclusively to estimate the mean time of travel of the 
tide From one point to another and the mean range at a given point. For ex- 
ample, if we look at the location of Port aux Basques and St. John's we note 
that the local values of Mo are 44 cm 14° and 35 cm 315°. Twice the amp 1i- 
tude gives 88 cm and 70 cm and the phase difference is equivalent to 2h 
2 min in time. The observed mean range at Port aux Basques and St. John's 


is 94 and 77 cm respectively; the mean time difference between high or low 
water between the two stations is 1 h 57 min. 
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Fig. 7. Observed amplitude ratio and difference in the Greenwich phase lag 


between the pairs K2 and S92, P] and K]. It is presented as a vector of 
length r (amplitude ratio) and inclined at an angle » (difference of phase 
lag). These values are obtained at stations where the constituents are 
separable; they are used afterwards for inferring K2 from S?2 and P] from ky 
at stations where they cannot be separated. Since the observed values depend 
on the geographical location, the measured vectors are presented in groups 
for the west coast and east coast. The tendency of the vector to 

bunch around the common value indicates that a mean value exists for the 
vector and the vector mean is presented in each case. The actually observed 
r, @ values are shown as rays in the larger diagram. Their vector average 
is shown in the smaller diagram as a single ray whose amplitude (r) and 
orientation (¢ , in degrees) are written out explicitly as "mean amplitude" 
and "mean angle". 
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6.1 Charts for Canada 


The charts (Fig. 8) show the tidal régime in all ocean waters adjacent 
to Canada and illustrate the effective response of these waters to the tide 
in the neighbouring Atlantic, Pacific and Arctic Oceans and to the local 
semidiurnal and diurnal tidal forces. In each the solid contours denote co- 
phase lines along which the vertical displacement associated with the spec- 
ific constituent is simultaneous. The phases are measured in degrees and 
correspond to the Greenwich phase lag in the Greenwich time zone (Z = 0). 
The tidal wave progresses in the direction of increasing phase. The time 
interval At corresponding to a given phase difference Ad may be obtained 
with the help of the following rule: 


Ag 
At = — hours (33} 
Oo 


where o is the frequency of the constituent given previously. For example, 
a phase difference of 20° between two contours is equivalent to: 


20 
At = ——  h = 41 min for Mo 
28.98 
20 
———-f) = OU Min or lon 2Z0umin for kK 
15.04 


Grossly one may say that 1° of phase is 2 min of time in the semidiurnal band 
and 4 min of time in the diurnal. The dashed contours are lines of equal 

amp 1itude and the units are centimeters. The contour interval is not regular 
and varies with the density of observations, the speed of change and the 
Space available to draw them. Therefore while studying the maps one must 
always watch for a sudden change in gradient, either in phase or in amplitude. 
It did not seem reasonable to draw regular and smooth contours in areas where 
there is virtually no information. On the other hand we had no fear of ex- 
trapolating contours in areas of more but still little information or of 
conjuring amphidromic structures wherever the Sparse data suggest one. 


We recall that the charts are of the vertical tide only; this represents 
an Incomplete picture of the tidal movement as it ignores the horizontal 
displacement or tidal current. In theory we can draw charts of Z or u with 
equal ease, whereas in practice the prospect as to when we will be able to 
draw charts of the tidal currents as actually observed is still remote. 


6.1.1 Nodes and antinodes 


We recall that the wavelength of K] iS approximately twice 


a 


Fig, 8. Cotidal chats for the. constituents Mol4So./ No, Ky “and’0)] in waters 
adjacent to Canada. A Miller projection is used for the geographic back- 
ground. The solid lines are contours of constant Greenwich phase lag; the 
dashed lines are contours of constant amplitude and the units degrees and 
centimeters. The contours are marked as 120° or 30 to indicate a Greenwich 
phase lag of 120° or a contour of 30 cm. The time zone used for all the maps 
is Greenwich (Z = 0). The spacing of the contours is not necessarily regu- 
lar. Question marks indicate areas where it was not possible to conjecture 
the type of motion present either because of a shortage of data or of their 
poor quality. The signs + or - help visualize the direction of increasing 
or decreasing values in areas where this is not evident. (a) My (b) So 
(c) No (d) Ky (e) Oj. 
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that of M2 in a sea of constant depth, so that an amphidromy in Mo (Fig. 8a) 
should imply an amphidromy in Kj} (Fig. 8d) further away from the reflecting 
wall. We find a well-defined amphidromy of Mo in Baffin Bay, off Clyde 
Inlet. The corresponding node of Kj falls on the Cumberland Peninsula and 
is degenerate. We note that friction displaces both nodes to the right of 
the reflecting wall, which in this case is the coast of Greenland at lati- 
tude 76°N in Nares Strait. The second amphidromic system for Mo is found in 
the Gulf of St. Lawrence and the corresponding node of Kj falls outside the 
gulf, south of Cape Breton Island and west of Sable Island. We find a second 
mininode of Mo in the northern portion of Northumberland Strait. 


There is evidence of amphidromic systems in Hudson Bay; they appear quite 
complicated because of the shape and shallowness of the bay. There are 
probably four semidiurnal nodes in a north-south direction; one lies at the 
latitude of Mansel Island, probably quite mangled because of the local top- 
Ography, one to the east of the Ottawa Islands, one south of Inoucjouac and 
finally one inside James Bay at the latitude of Moar Bay. Friction has 
pushed them well to the right of the reflecting wall (the southern portion 
of Hudson Bay and James Bay) and the node in James Bay is completely degen- 
erate. Because of the width of Hudson Bay an additional transverse node 
exists east of Churchill. The diurnal counterpart of the system described 
consists of a node south of Coats Island and another (degenerate) at Cape 
Jones. 


It is almost certain that an amphidromy exists in the area off Coppermine 
with another likely in Foxe Basin where shallow depths imply short wavelengths. 
There probably exists an amphidromy west of Sachs Harbour. We see no evidence 
of nodes in the North Pacific off British Columbia. 


Antinodes are never as obvious as nodes in cotidal charts, nevertheless 
between two nodes there should always be an antinode. A spectacular example 
is found in Hudson Strait where it is marked by extreme amplitudes in the 
semidiurnal band because of resonance and a near absence of cophase lines 
indicating simultaneous tides over a large area. The corresponding diurnal 
antinode is found in the Labrador Sea. A more emphatic diurnal antinode is 
to be found in Committee Bay west of Fury and Hecla Strait where the diurnal 
tides are as large as those on the west coast in Georgia Strait. 


Because of the numerous openings of Canadian basins to the ocean, we 
encounter areas where the partial waves from various directions come together, 
areas where "tide meets tide". With the help of cotidal charts we may be 
quite precise when we speak of "tide meets tide". It is the area where the 
phase of two waves coming from different directions is the same; in a zone 
little affected by friction it would also mark the point where the amplitude 
is maximum. We note for example that tide meets tide to the north east of 
Queen Charlotte Islands for the case of Mo, but that the maximum amplitude is 
further south (Fig. 9). The actually observed tide is made up of the super- 
position of all the harmonic components and the zone of "tide meets tide" in 
the observed situation may move from day to day and from neap to spring tide. 
The phenomenon can be observed in many areas including: Kane Basin in Nares 
Strait for the semidiurnal tide (the diurnal tide progresses from the Atlan- 
tic northward to the Arctic Ocean), the Strait of Belle Isle for both tides 
and in Hecate Strait. 
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The largest but least obvious example of the phenomenon is found west 
of the arctic archipelago; the archipelago acts as a complicated conduit 
through which the Atlantic tide makes its way westward to meet the tide in 
the Arctic Ocean. The opposite occurs inside the shallow Coronation Gulf 
which is penetrated by the arctic tide as far as King William Island. We 
shall see more examples when we scrutinize the larger scale maps. 


6.1.2 Areas in need of further exploration 


We have indicated in the map of M2 (Fig. 8a) those regions 
which are essentially devoid of tidal observations; naturally they fall in 
the arctic, an immense and quite inaccessible area which will be known only 
with the help of a systematic and rational program of exploration. A list 
of the regions which need to be explored are: the Canadian shore of the 
Arctic Ocean from Banks Island to Greenland; inside the arctic archipelago; 
Amundsen Gulf; Coronation Gulf, McClintock Channel and the Gulf of Boothia; 
Foxe Basin and the northern entrances to Hudson Bay; the west coast of James 
Bay from Cape Henrietta Maria to Sand Head; the east coast of Ungava Bay; 
the Bay of Fundy and the east coast of Hudson Bay. I mention this region 
last because it presents some special difficulties. Study of the records of 
Inoucjouac indicates that the locally weak and fickle tides are very much 
affected by the formation of an ice cover; the range varies by 60% and the 
mean time of arrival of the tide may vary by one hour from season to season. 
Therefore some special care must be taken before developing a plan for mea- 
Surements there. 


6.2 Patterns of tides 


We can draw contours of the amplitude and phase of the constituents o, 
as well as of patterns of tide (Fig. 9). Such a map shows at a glance the 
various types of tide present and where they are to be found. The nodes 
delineate zones of rapid changes in patterns because of the steep gradients 
in amplitude they imply; therefore we shall meet pure D or pure S tides at 
the sites of the nodes we noted previously. We see such areas in the shallows 
around the Magdalen Islands, in Northumberland Strait, south of Cape Breton 
Island, in Hudson Bay and in Baffin Bay. 


The tide is strongly semidiurnal south of the Maritimes and Newfoundland, 
in Hudson Strait and in the Bay of Fundy. It is mainly semidiurnal in the 
Labrador Sea, in the northern section of the Gulf of St. Lawrence and in the 
arctic archipelago. It is strongly diurnal in Committee Bay. It has a 
weakly mixed character along the Pacific west coast, which indicates the ex- 
istence of an appreciable diurnal inequality. It acquires a more mixed 
character in Juan de Fuca and Georgia Straits because of a degenerate semi- 
diurnal node around Victoria. 


In all areas of the map where the character of the pattern is left of 
.6M, we may compute lunitidal intervals and ranges unambiguously. 
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7 LARGE SCALE MAPS 


7.1 British Columbia 


The small scale maps (Fig. 8) indicated that the tidal motion in our 
section of the Pacific is relatively simple and consists of a counterclock- 
wise progression of the crest of high water with increased amplitude north- 
ward. The large scale (Fig. 10) shows additional details east of the Queen 
Charlotte Islands and east of Vancouver Island. The map for each constituent 
is given in two sheets; one covering the northern section from Alaska to 
Vancouver Island, the other the southern section around Vancouver Island. 
The tidal motion gets more complicated than in the open Pacific and consists 
essentially of the breaking up of the tidal wave into northern and southern 
branches. A finer scale motion consists of the progress of the tide up the 
numerous fjords indenting the coast in which the details of the motion are 
of local interest only. 


In the north the semidiurnal Pacific wave enters Hecate Strait from Dixon 
Entrance and from the broader Queen Charlotte Sound. There is a zone of 
simultaneous phase between Rose Point and Skidegate Inlet with a more defined 
zone of maximum amplitude off Skidegate Inlet. This holds as well for So and 
No. This area defines a zone of "tide meets tide" where large tides are 
usually accompanied by weak currents. The progress of the tide inside the 
various inlets is indicated wherever information is available. 


In the south the local geography creates further complication in the 
progress and development of the tide. Once again the Pacific semidiurnal 
wave divides into two branches, the northern one entering through Queen Char- 
lotte Strait and the southern one through Juan de Fuca. Because of the 
shallowness of Juan de Fuca a semidiurnal node is formed at the longitude of 
Victoria-Esquimalt which is degenerate. There is a delay of 23 hours between 
the entrance of the tide in Juan de Fuca and its passage by Victoria. The 
degenerate amphidromy retards it further and then also the Gulf Islands, so 
it emerges into Georgia Strait five hours after its entrance into Juan de 
Fuca. There is some complicated motion inside the Gulf Island caused by the 
local geography, which has been described elsewhere (Godin 1976) but in 
Georgia Strait itself the tide is effectively quasi-simultaneous as far as 
Campbell River. Here we see a good example of a zone of antinode, the con- 
stricted geography in the north acting as an effective barrier to the semi- 
diurnal wave entering through Juan de Fuca. The northern branch of the semi- 
diurnal wave progresses normally as far as Johnstone Strait. Further south 
it is effectively blocked by the system of narrows north of Georgia Strait. 
We note extreme delays of the order of two hours within a span of 10 nautical 
miles in Seymour Narrows and Hoskyn Channel. It enters Georgia Strait with 
a phase equal to that of the wave coming from the south indicating the limit 
of progress of the northern wave. 


_ The Mg wave in the Pacific has an amplitude of 90 cm in the southern sec- 
tions of British Columbia increasing to 110 cm near Alaska. The Mo signal is 
further amplified in Hecate Strait and reaches a maximum of 200 cm off 
Skidegate. In the south the reverse occurs in Juan de Fuca where the Mo 
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Fig. 10. Cotidal charts for the constituents M2, So, No, Ky] and 07 in 
waters adjacent to the northern and southern portions of British Columbia. 
(a) Mo (b) So (c) No (d) Ky (e) 0}. 
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amplitude decreases from 90 cm at the mouth to 35 cm off Victoria-Esquimalt. 
The degenerate amphidromy creates an east-west gradient in amplitude towards 
Admiralty Inlet. Thereafter the amplitude increases rapidly through the 

Gulf Islands and the transverse gradient disappears eventually in Georgia 

Strait. In the northern entrance the Mo amplitude increases regularly from 
100 to 130 cm in the entrance of Johnstone Strait. Losses diminish it after 
and it gets to 95 cm in the northern entrance to Georgia Strait. A similar 


description applies equally to the other semidiurnal waves S2 and No Lig, 
Ob and -c)- 


The modification of the diurnal Pacific tide by the Queen Charlotte Is- 
lands and Vancouver Island is quite similar to that for the semidiurnal tide; 
a difference is that the diurnal node equivalent to the node near Victoria- 
Esquimalt lies outside Juan de Fuca. The Kj amplitude has a zone of amp 1 i- 
tude larger than 50 cm in the northern portion of Hecate Strait. Over the 
Pacific all along the Canadian coast it has an amplitude of the order of 
40 cm. It increases to 50 cm in Queen Charlotte Strait. In Juan de Fuca it 
increases steadily from 45 cm to 80 cm in its entrance to Georgia Straits 
There is no local minimum off Victoria-Esquimalt because of the absence of a 
diurnal node. Amplitudes larger than 90 cm are found around Texada Island. 
There is indication that the southern branch of the diurnal wave diffuses 
somewhat through the narrows north of Georgia Strait. 


The patterns (Fig. 11) indicate that the tide in the Pacific has a slight- 
ly mixed character. It becomes more normal in Hecate Strait and although the 
number of contours is quite large, the tide over the whole area has a rather 
uniform character. This changes abruptly in Juan de Fuca because of the 
semidiurnal node there and it has a character higher than .6M (mixed) from 
Port Renfrew all the way to Georgia Strait and beyond Seymour Narrows. The 
calculations of conventional lunitidal intervals in that region cannot be 
too successful because of the presence of intervals during which the tide is 
diurnal. 


7.2 Gulf of St. Lawrence, Cape Breton and Newfoundland 


The semidiurnal tide (Fig. 12a, Db, c) is dominated by the amphi dromy 
west of the Magdalen Islands. The node is strongly distorted by the presence 
of the islands and the very few data available seem to indicate that east of 
the islands there is still a marked semidiurnal signal, thus inducing abrupt 
gradients in the amplitude contours. We find a second node at the northern 
end of Northumberland Strait. It may be the same node encountered in the 
Magdalen Islands which has been split by Prince Edward Island or it may have 
been created by the local bathymetry. It is still intact for Mp and No in 
spite of the narrowness and shallowness of the channel] but appears to have 
been pushed against the coast in the case of S2 (Fig. 120) 


The movement of the Mp tide in the Gulf of St. Lawrence, because of eee 
presence of the Magdalen node, consists of the counterclockwise motion of ai 
crest of high water which enters Cabot Strait then progresses along the = 
coast of Newfoundland, the Québec northshore, Gaspé Peninsula and the north- 
eastern coast of Prince Edward Island. The motion inside Northumberland 
Strait consists of a wave entering from the south progressing up the strait 
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Fig. Il. Patterns of the tide in waters adjacent to the northern and 
southern portions of British Columbia. (a) Northern portion (b) Southern 
portion. 
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with increasing amplitude until it meets the region of the node, where there 
is a rapid decrease in amplitude and marked delays in the progress of the 
tide. 


The semidiurnal tide is nearly simultaneous along the east coast of 
Newfoundland while along its southern shore, it starts progressing towards 
Cabot Strait taking about two hours to reach it. Past the Gaspé Peninsula it 
takes a progressive character and travels towards Québec with continuous am- 
plification. The maximum amplitude is found between Isle aux Couldres and 
the eastern tip of Ile d'Orléans. There the channel splits and the depth 
decreases, friction takes over and the tide acquires an increasing shallow- 
water character. It takes eleven hours to progress from Québec to Sorel at 
the western extremity of Lake St. Peter beyond which it is no longer detect- 
able. 


The Kj tide (Fig. 12d) inside the Gulf of St. Lawrence presents a simpler 
picture, because its node lies outside and it treats the Gulf area as far as 
Pointe des Monts and the Strait of Belle Isle as an integral part of the 
Atlantic Ocean. There is a time delay of over 5 hours between the diurnal 
peak at Port aux Basques and North Sydney. It meets the Belle Isle branch 
of the wave well inside the gulf, just north of Anticosti, because the Strait 
of Belle Isle is also quite permeable to the diurnal signal. On the other 
Side it progresses towards the Magdalen shallows and splits in two branches 
north and south of Northumberland Strait. The diurnal tide is virtually 
Simultaneous along the east coast of Newfoundland (Fig. 12e) and like the 
semidiurnal tide it starts progressing along its southern shore towards Cabot 
Strait taking nearly 7 hours. Because it has a small amplitude in the open 
Atlantic and in the gulf area, the diurnal tide is of secondary importance. 
It is amplified in Northumberland Strait and off the New Brunswick coast, 
creating local diurnal inequalities. The diurnal tide takes a progressive 
character in the St. Lawrence estuary increasing in amplitude as far as the 
Ile d'Orléans. Thereafter it is increasingly delayed and distorted by fric- 
tion; it is no longer detectable past Louiseville in the middle of Lake St. 
Pierre having taken 11 hours to reach it from Québec (Fig. 13). 


Semimonthly tides are perceptible at Montréal (Fig. 14) although there 
are no regular tides as we know them near the seashore. This tide is created 
mainly by the nonlinear interaction of Mz and So and reflects the occurrence 
of spring and neap tides downstream by rises and falls in the local level. 
The effect becomes apparent at Cap a la Roche and reaches a maximum at Trois- 
Riviéres where it exceeds in amplitude the regular tide. It is quite faint 
but still perceptible in Montréal where effects of the river discharge are 
also apparent. 


The patterns in the geographical area which we investigate are quite com- 
plicated because of the nodes (Fig. 15) except along the east coast of New- 
foundland where they are almost uniformly of the N category. South of New- 
foundland they move rapidly from the .4S (strong semidiurnal) category to the 
N (normat) category in Cabot Strait. The patterns in the northeastern sec- 
tion of the gulf hover around the boundary of the N and M category. From 
Anticosti to the Magdalens they shift rapidly from .2M (slight diurnal in- 
equality) to 2.0D (diurnal). They are mixed along the north shore of Prince 
Edward Island and they pass through the M and D category inside Northumber- 
land Strait. In the estuary the tide reverts to slightly mixed then to the 
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Fid. (3. Cotidal charts for the constituent Mo,.S5, No, Ky and Oy in’ the 
upper St. Lawrence River. (a) Mo (b) Sa (c) A> 4) Ky S OT. 
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14. Creation of a semimonthly tide in the upper reaches 


of the St. Lawrence River caused by the succession of spring 
and neap tides downstream (from Godin 1979). 
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Fig. 15. Tidal pattern around the Gulf of St. Lawrence and the upper 
estuary of the St. Lawrence River. (a) Gulf of St. Lawrence (b) Upper 
estuary. 
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normal category to finally assume a semidiurnal character past Pointe des 
Monts. It is semidiurnal inside the St. Lawrence River, but as the semi- 
diurnal signal gets damped more rapidly than the diurnal it gradually ac- 
quires a stronger diurnal inequality. 


7.3 Bay of Fundy 


Although this area is world famous the actual data we possess about 
the tide over the body of the Bay of Fundy is quite inadequate. This re- 
sults exactly from the fact that large tides imply foreshores that dry for 
miles and ranges which tax most instruments. The charts which we present 
(Fig. 16) should be considered as schematic at best. The wave from the Gulf 
of Maine divides and its eastern branch enters the Bay of Fundy. There 
seems to be some diffraction around Grand Manan Island. 


Inspection of the charts for Mz, So and No (Fig. 16a, b, c) indicates 
that No is the true resonating semidiurnal constituent in the Bay of Fundy. 
The No tide is almost simultaneous over the whole bay and its amplitude in- 
creases threefold between Grand Manan and Cape d'Or. In contrast M2 only 
doubles over the same stretch and So less than doubles; Mp takes 1/2 hour to 
reach Cape d'Or and Sp about an hour. The largest amplitudes are found in- 
Side Minas and Cumberland Basin in which submerged gauges were deployed, thus 
increasing the quality of the information. The tide takes nearly an hour to 
penetrate Minas Basin in which it is further amplified. This time Mo seems 
more favored because it increases a further 50% while No gains a mere 10%. 


The diurnal tide (Fig. 16d, e) is similar in profile to the semidiurnal 
tide except that it falls far away from resonance and is little amplified in 
the upper regions of the bay. 


8. USE OF COTIDAL CHARTS FOR PREDICTION 


The contributions of the five constituents illustrated in the charts, in 
conjunction with P] and Kg, usually amount to over 90% of the tide observed; 
therefore they can be used to prepare predictions of the tide whenever a 
high degree of precision is not necessary, or at points in the ocean where 
the tide has never been observed. The cotidal charts represent an intelli- 
gent extrapolation of the information available and although they may be 
deficient in some details, they constitute the best guess as to what is 
happening. 


Should a prediction be required at an ungauged site, its position is 
located on the charts and the corresponding values of M2, So, No, Kj and 07 
are obtained by interpolation between the cotidal lines. The additional 
values of Py and Ko are inferred by using the regional values of r and 4% and 
applying them to Kj and So. Then it must be decided in which time zone the 
prediction is to be made since the charts are in Greenwich time. In a sci- 
entific application it is preferable to retain the Greenwich time because 
the switch to other time zones is an unending source of error and confusion. 
If there is reason to prefer to work in some other time zone, the Greenwich 
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phase lags of the constituents have to be first transformed to the chosen 
time zone with the help of transformation (4). At this stage we are provid- 
ed with a set of amplitude and phase (A, g) for the seven constituents men- 
tioned. This is not sufficient. We also need the astronomical argument V 
of the constituent, which is derived from astronomical knowledge and the 
nodal corrections (f, u) derived from tidal theory, which modify the ampli- 
tude A and the argument V of each constituent. In addition, although we 
always write a cosine in (1) for the contribution of all constituents, theory 
informs us that diurnal constituents actually enter as sines, while semi- 
diurnal constituents contribute as cosines. In addition their contribution 
may be positive or negative, and this additional information is also provided 
by theory. 


The calculation of angles, when there is a possibility that we will be 
consistently faced with angles containing very large multiples of 360°, is 
usually done with the help of a special angular unit called "cycles". An 
angle is measured by the portion of a complete circle it covers; therefore 
one cycle corresponds to a full circle or 360° in degree measure or 2t 
radians in radian measure. Thus 45° in cycles becomes: 


45/360 = 1/8 = .125 cycle 


The advantage of cycle measure is that if we had to multiply 45° by 75, we 
obtain 3375° which has to be reduced to 135°. The same calculation in 
cycles becomes: 


io Xe. 125 = 0. 32s cvele —=.,375. cycle (1357) 


The reduction is done simply by dropping any whole number of cycles in the 
final result since whole cycles represent complete turns around the circle. 
Henceforth we shall use cycles in the elaboration of formulas for predictions 
with the understanding that the argument of the sine or cosine is expressed 
in degrees or radians in the ultimate step. 


As an example we set up a scheme of predictions for Churchill starting 
at 00 hour (Central Standard Time) June 5, 1980. In order to do this we 
need two additional tables: one giving the number of days elapsed since the 
beginning of the year (Table 4), the other the astronomical argument V of 
each constituent at 00 hours GMT at the beginning of the year and the mean 
nodal corrections (f, u) for that (Table 5) year. The latter tables were 
prepared by the German Hydrographic Institute and cover the years 1900-1999 
(Anon. 1967). Table 4 indicates that June 5 is the 157th day of 1980, a leap 
year. Table 5 gives the argument V for January 1 at 00 hour; we therefore 
have to increase it by adding the number of hours elapsed between January 1 
and June 5 multiplied by the angular speed of the constituent in order to 
adjust to the initial time of the prediction. The table actually gives 
Vo + u, i.e. automatically incorporating the nodal correction, so we simply 
add the adjustment to V9 + u to get V + u at the desired date. Then one 
multiplies the amplitude A by f in order to obtain the adjusted amplitude of 
the constituent during the interval of predictions. 


8] 


Table 4. Number of days elapsed since 
the beginning of the year (from Godin 
W792. 0. 16). 


Number of days elapsed 


Number in the year at the 

mm Month of days beginning of the month 
01 January 31 000 

02 February 28 (29) 031 

03 March 31 059 (060) 

04 April 30 090 (091) 

05 May 31 120 (121) 

06 June 30 151 (152) 

07 July 31 181 (182) 

08 August 31 212 (213) 

09 September 30 243 (244) 

10 October 31 273 (274) 

Val November 30 304 (305) 

12 December 31 334 (335) 


Table 6 shows this task of preparation for Churchill. The Zg in the formula, 
e.g. (1) and Table 6, stands for the reference level. It is usually unknown 
in unexplored areas and is then replaced by zero, so that the predictions 

are about mean sea level. If one wants to obtain the hourly predicted levels, 
t is increased by increments of 1 hour. If the time and height of high and 
low water are desired, one looks for the zeros of the derivative of z(t); a 
practical scheme for such a search could be the one Suggested by Godin and 
Taylor (1973). The map of patterns indicates if the local tide will be of 
semidiurnal, mixed or diurnal character. | 
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Table 5. Tables of f and Vy + u between 1900 and 1999 for 
Ose PTs Kys Nas Mo. Ss) and Ko from tables prepared by the 


German Hydrographic Institute (Anon. 1967). 
te u Migs U £ U 
1900 0.950 -010.88 Ef ie 1950 15183 +000.04 
071 0.893 -009.99 OSV! 51 Moss -002.48 
02 0.847 -006.99 1O4ue 52 1.154 -004.91 
03 0.815 -002.50 296.6 53 eid -007.67 
04 0.806 +002.19 04253 O,-= 145,555 54 1.068 -009.94 
05 0.830 +006.07 ete 55 On EOOse2 
06 0.884 +008.91 225.6 56 0.948 -010.43 
O7 0.949 +i) Oven a 528.4 57 0.880 -009.11 
08 1.009 +011.08 069.8 58 0.832 -006.54 
09 061 +009.74 144.1 59 0.812 -002.38 
1910 eS +007.43 242.8 1960 0.816 4002.72 
1 1.160 +005.17 341.5 61 0.842 +007.15 
12 4c +002.87 080.2 62 0.887 +009.72 
1) WS AS -000.11 M52. 63 0.951 TOMO o> 
14 1.169 -003.41 250.5 64 #2022 +010.41 
115) es ios) -006.03 4348.9 65 ee +009.61 
16 Absaliailiezs -007.93 (Osta! 50) 66 1s WeO +007.65 
cs 1.052 -009.69 161.8 67 ee +004.71 
18 0.990 2014.03 261.5 68 4, 184 +001.83 
19 0.933 -010.86 002.7 69 1.188 -000.61 
1920 0.877 -008.69 105.8 1970 167 -003.28 
21 0.828 -005.09 185.0 74 1.136 -~006.34 
Bo 0.803 -000.87 290.3 Jo i, 162 -008.82 
23 0.813 +003.40 035.5 73 1.053 =010.10 
24 0.854 +007.30 140.4 74 0.985 -010.55 
25 0.910 +010.15 218.9 75 0.915 -010.32 
26 0.967 4014.21 320.9 76 0.862 -008.66 
27 1.026 +010.40 061.1 17 0.827 -004.98 
28 OSS +008 .66 160.4 78 Oren -000.01 
29 139 +006.83 234.2 79 0.818 +004.72 
1930 165 +004.61 332.9 1980 0.854 +008.03 
31 em +001.54 070.9 81 0.914 +009.91 
32 178 SOO. te 168.6 82 0.984 +010.079 
Se ees -004.31 241.6 83 1.044 +010.57 
34 le 156 -006.48 340.5 84 1.090 +008.86 
35 1.084 -008.79 079.1 85 1.135 +006.19 
36 1.030 01070 178.2 86 A473 +003.63 
57 0.975 S011 510 253.4 87 1.186 +001.26 
38 0.911 -009.83 355.7 88 rc MP 001 62 
39 0.850 -007.34 099.2 89 are -004.86 
1940 Os8i2 -003.88 2 OS 1990 leeulhore -007.43 
41 0.807 +000 ail 283.6 91 1.086 -009.05 
42 0.831 +005 .20 029.3 92 Teor -010.26 
43 0.874 +008.97 134.1 93 0.954 -010.95 
44 0.926 +010.78 2509 94 0.898 -010.11 
45 0.988 +010.67 312.4 95 0.850 007.17 
46 1.057 +009.67 052.4 96 0.815 =002.76 
AF tee es +008.34 Tee 97 0.805 +001.85 
48 sri AS) +006.22 ZO st 98 0.828 +005 77 
49 1.165 +003.15 323.5 99 0.881 +008.77 


. «ae © -«° “Ss Gabe key ay. tery .i6) 0: 
= 6) oS) SES ko 1S Oo 


CU ee en aE aE cls Wr GS 
Nw fF FW UA Oo WW PP 


NM) Srey WA Gy Wy Vi OS 


of el. eee sats eC te 
= Ora © OF © OM 
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SN! GA JS ROY. OP sr ie 
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—> WW 


1900 


1910 


1920 


1940 


U 


+000. 66 
+000.48 
+000.34 
+000.22 
-000.03 
-000.36 
-000.56 
-000.58 
-000.56 
-000.59 


-000.56 
-000.35 
-000.08 
+000.08 
+000.16 
+t OOOKmm 
+000.58 
+000.71 
+000.65 
+000.54 


+000.49 
+000.37 
+000.09 
-000.24 
-000.44 
-000.49 
-000.56 
-000.66 
-000.63 
-000.42 


-000.19 
-000.07 
+000.03 
+000.26 
+000.52 
+000.63 
+000.61 
+000.59 
+000.61 
+000.50 


+000.21 
-000.11 
-000.28 
-000.39 
-000.55 
-000.70 
-000.66 
-000.47 
-000.30 
-000.22 
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Habbe 5. .2(Gontd. ) 
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350.5 51 
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549...7 1960 
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S5alit 1970 
350.3 71 
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Table 5. (Contd.) 


f u Vig £ u Vert 

1900 ~—-0.973 +008 .66 018.8 1950 es -000.33 009.8 
01 0.935 +007. 46 017.4 51 1.109 +002.24 1241 
02 0.904 +005.19 014.9 52 1.096 +004.64 014.2 
03 0.885 +002.07 011.5 53 1.073 +006.67 017.0 
04 0.884 -001.42 007.8 K, = 165.555 54 1.043 +008.14 018.3 
05 0.899 -004.65 005.3 55 1.006 +008.83 018.7 
06 0.928 007.1" 002.6 56 0.967 +008.54 018.2 
07 0.966 -008.52 001.0 57 0.929 +007.16 017.5 
08 1.005 -008.83 000.4 58 0.900 +004.73 014.9 
09 1.042 -008.17 001.8 59 0.884 +001.51 011.4 
1910 1.073 -006.72 003.1 1960 0.885 -001.97 007.7 
11 1.095 -004.70 004.8 61 0.903 -005.11 005.3 
1 1.108 ~002754 007.0 62 0.934 -007.41 002.8 
13 tonae +000. 26 010.3 63 0.972 -008.64 001.3 
14 1.106 +002.80 012.6 64 Ot -008.79 000.9 
15 1.091 +005.13 014.7 65 1.047 -007.98 002.5 
16 1.067 +007.06 016.4 66 POH -006.43 003.8 
17 1.035 +008.37 018.4 67 1.098 -004.33 005.6 
18 0.997 +008 .86 018.7 ; 68 rea tace -001.90 007.8 
19 0.958 +008.33 017.9 69 ae +000.67 Cat 
1920 0.922 +006.70 016.1 1970 1.105 +003.20 013.4 
21 0.895 +004.05 014.2 71 1.088 +005.48 015.5 
22 0.883 +000.73 010.6 72 1.063 +007.32 017.1 
23 0.888 =002273 006.9 73 1.029 +008.51 019.0 
24 0.909 -005.72 003.7 74 0.991 +008.84 019.1 
25 0.942 =007.79 002.3 75 0.952 +008.14 018.2 
26 0.981 -008.77 001.1 76 0.916 +006.34 016.1 
27 1.020 -008.68 001.0 717 0.892 +003 .55 014.1 
28 1.055 -007.69 001.7 78 0.882 +000. 16 010.5 
29 1.082 -005.99 004.2 79 0.890 -003.26 006.8 
1930 Feo4 -003.81 006.1 1980 0.914 -006.12 003.7 
31 as -001.33 008.4 81 0.948 -008.02 002.5 
a tA $001.25 010.7 82 0.987 -008.82 001.5 
33 $i102 +003.74 013.9 83 1.026 -008.58 001.5 
34 1.083 +005.94 015.9 84 1.060 -007.46 002.4 
35 1.056 +007.65 017.4 85 1.086 -005.67 004.9 
36 1.021 +008.67 018.1 86 1.104 -003.42 006.9 
37 0.982 +008.78 019.0 87 ee -000.91 009.2 
38 0.943 +007 .83 017.8 88 1.110 +001.67 011.5 
39 0.910 +005.79 015.5 89 1.099 +004.12 014.7 
1940 0.888 +002.82 012.3 1990 1.079 +006 .26 016.6 
41 0.882 -000.63 009.6 1 1.050 +007.87 018.0 
42 0.894 -003.97 006.0 92 1.015 +008.75 018.7 
43 0.921 -006.64 003.1 93 0.976 +008.70 019.4 
44 0.957 -008.30 001.2 94 0.937 +007.57 018.0 
45 0.996 -008.85 001.4 95 0.905 +005.36 015.5 
46 1.034 -008.40 001.7 96 0.886 +002.28 012.2 
47 £066 =007 201 002.7 97 0.883 -001.20 009.5 
48 1.091 -005.19 004.4 98 0.898 -004.46 006.0 


49 TOG -002.87 007.5 99 0.926 -006.98 003.62 
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Tapre’ 5. -<( Corned.) 


i u Viot 1 a u Vo+u 

pe Rees +001.83 065-5 1950 0.964 +000.12 216.0 
01 1.021 +001.47 OTT 2 51 0.968 +000.55 228.5 
02 1031 +001.20 089.0 52 0.968 +001.20 B51. 
03 1.038 +000.34 100.1 53 0.980 +001.92 216.4 
04 gO: -000.33 Tt SS N, = 245.655 54 0.993 +001.93 228.5 
05 1.033 -000.73 O85 aif 55 1.000 +002.02 240.6 
06 1.029 -001.55 096.9 56 1.015 +002.10 25207 
07 1.014 -001.99 108.5 57 1.028 +001.46 226.7 
08 1.004 “060697 120.5 58 1.030 +000.93 238.2 
09 0.994 062/20 094.9 59 T5057 +000.50 249.8 
1910 0.977 -001.93 107.2 1960 1.040 ~000.48 260.8 
11 0.970 -001.15 120.0 61 1.030 -001.14 avin 
12 0.967 -000.73 132.4 62 1.025 -001.46 246.5 
13 0.960 +000.03 107.8 63 1.016 -002.09 257.9 
14 0.966 +001.02 120.8 64 0.998 -002 522 269.8 
15 0.975 +001.35 M3562 65 0.989 -001.86 244.7 
16 0.980 +001.79 145.6 66 0.980 -001.79 256.8 
17 0.995 +002.24 120.7 67 0.967 SOO1025 269.4 
18 1,009 +001.95 13204 68 0.965 -000.34 282.3 
19 1.015 +001.73 144.2 69 0.966 +000.09 257.3 
1920 tROaF +001.56 156.1 1970 0.964 +000.80 270.1 
a 1.036 +000.75 129.9 71 0.974 +001.64 282.9 
22 1.034 +000.11 141.2 72 0.986 +001.78 295.1 
23 1.036 “O06034 152.8 ie 0.992 +002.04 269.9 
24 1,034 “000023 164.0 74 1,008 +002.25 282 2 
25 1.020 —OONCT? 138.0 75 1.022 +001.73 293.7 
26 Neots -001.88 149.9 76 1.026 +001.29 305.3 
mt 1.002 -002.26 161.6 77 1.036 +000.90 279.5 
28 0.984 2602012 ‘73.8 78 1.040 -000.07 290.5 
29 0.975 -001.46 149.0 79 P0353 -000.75 301.9 
1930 0.970 SOO1G4 161.4 1980 1.030 -001.14 313.5 
31 0.961 -000.42 174.1 81 i025 -001.88 287.4 
22 0.964 +000.58 87 «1 82 1.006 002213 299.2 
53 0.971 +000.98 162.1 83 0.996 -001.93 311.4 
34 0.974 +001.53 WT4AS7 B4 0.987 -002.01 505.5 
35 0.987 +002.12 1873 85 0.972 -001.59 298.4 
36 1,001 +001.98 199.2 86 0.968 -000.77 311.2 
Sif 1.008 +001.92 Siok 87 0.966 -000.37 323.6 
3 13027 +001.86 135 27 88 0.962 +000.37 336.4 
39 1,032 4001.13 497 20 89 0.969 +001.29 S169 
#900 95033 +000.53 208.4 1990 0.9779 +001656 aah 
41 1,038 +000.08 182.6 91 0.985 +001.96 aro 
42 1.038 -000.87 193.7 92 1.000 +002.32 349-0 
43 1.026 -001.48 205.1 93 +15 +001.92 323.2 
441.019 001.71 216.9 ga aoe et Bees? 
45 1.009 -002.22 191.0 95 1.033 deli  sitahiy 
46 04991 -002.22 203.0 96 1.039 +000.34 hid 
ie 0.982 -001.71 215.5 97 15.035 -000.34 351.6 
48 0.974 -001.50 227.8 og 1.034 -000.78 es 
49 0.963 -000.85 203.0 99 1.028 -001.60 354.4 


1900 


1910 


1920 


1930 


1940 


Uu 


+001.99 
+001.65 
+001.09 
+000.40 
-000.30 
~000.95 
2001351 
-001.93 
~002.15 
-002.07 


-001.74 
-001.27 
-000.68 
+000.05 
+000.81 
+001.42 
+001.82 
+002.07 
+002.12 
+001.89 


+001. 43 
+000.81 
+000.14 
-~000.54 
-001.18 
-001.72 
-002.06 
-002.12 
-001.94 
-001.60 


-001.09 
-000.38 
+000.39 
+001.05 
+001.56 
OOo 
+0026 15 
+002.06 
aOLO II aft 
+001.20 


+000.58 
-000.11 
-000.81 
-001.44 
-001.88 
-002.08 
-002.07 
-001.87 
-001.45 
-000.79 


Fable: 5. _ Contd. ) 
Veer U 

008.3 1950 
MOSiant Bil 
208.9 52 
309.0 53 
O4920s My 255.555 54 
124.8 55 
uA 10) 56 
325.3 57 
065.8 58 
142.3 59 
DUNS) 5 isk 1960 
344.6 61 
085.9 62 
163.0 63 
264.5 64 
005.9 es 
107.1 66 
183.7 67 
284.5 68 
025.0 69 
125.3 1970 
PRO al al 
SOMA V2 
Ova 2 5) 
WA, 3 74 
277.2 75 
B76 76 
058.3 El 
ton Tus} 
235.9 79 
aT. 2 1980 
078.6 81 
eOne 82 
257 2 83 
Biot! 84 
099.6 85 
200.5 86 
ZOMS toh 87 
Onan 88 
es 89 
PV aS 1990 
293.3 91 
OFS 3 92 
133.5 93 
233.8 94 
309.9 95 
OSHONS 7/ 96 
ae) 97 
BARI 98 
329.9 99 
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u 


-000.05 
+000.62 
+001.23 
+001.76 
+002.09 
+0025 13 
+001 .92 
+001 .53 
+000.99 
+000. 33 


-000.40 
-001.09 
-001.62 
-001.95 
-002.10 
-002.06 
-001.74 
-001.16 
-000.49 
+OOOR 17 


+000.85 
+001.49 
+001.94 
OO) 
+002.05 
+001.81 
OOnreort 
+000.74 
+000.02 
-000.69 


-001.29 
-001.75 
=002.°05 
-002.15 
-001.95 
-001.49 
-000.92 
-000.29 
+000. 43 
+001.16 


+007 270 
+001.99 
TOO LO 
+002.01 
+001.68 
+001.12 
+000.45 
-000.25 
-000.91 
-001.48 
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Table 5. (Contd.) 


£ u Vor f u 

1900 0.999 -000.12 359.9 1950 1.002 +000.01 
01 0.999 -000.10 359.9 51 1.002 ~000.04 
02 0.998 -000.06 359.9 52 TeOOr -~000.07 
03 0.998 -000.03 000.0 53 1.001 ~000.10 
04 0.998 +000.02 000.0 S, = 273.555 54 1.001 -000.12 
05 0.998 +000.06 000.1 55 1.000 -000.13 
06 0.998 +000.09 000.1 56 0.999 -000.12 
07 0.999 +000.12 000.1 57 0.999 -000.09 
08 1.000 +000.13 000.1 58 0.998 -000.06 
09 1.001 +000.12 000.1 59 0.998 -000.02 
1910 1.001 +000.10 000.1 1960 0.998 +000.02 
qi 1.002 +000.08 000.1 61 0.998 +000.06 
12 1.002 +000.04 000.0 62 0.999 +000.10 
1) 1.002 000.00 000.0 63 0.999 +000.12 
14 1,002 -000.05 000.0 64 1.000 +000.13 
15 1,002 -000.08 359.9 65 1.001 +000.12 
16 1.001 -000.11 359.9 66 1.001 +000.10 
Az 1.000 =000,12 359.9 67 1.002 +000.07 
18 1.000 -000.13 359.9 68 1.002 +000.03 
19 0.999 000.11 359.9 69 1.002 -000.01 
1920 0.998 -000.09 359.9 1970 1.002 -000.05 
25 0.998 -000.05 359.9 71 V002 -000.09 
22 0.998 -000.01 000.0 72 1.001 -000.11 
23 0.998 +000.03 000.0 73 1.000 -000.13 
24 0.998 +000.07 000.1 74 1.000 -OCome 
25 0.999 +000.10 000.1 75 0.999 -000.11 
26 0.999 +000.12 000. } 76 0.998 -000.08 
27 1.000 +000. 13 000.1 77 0.998 -000.04 
28 1.001 +000.12 000.1 78 0.998 000.00 
29 1.001 +000.09 000.1 79 0.998 +000. 04 
JSEIO O02 +000.06 000.1 1980 0.998 +000.08 
31 1.002 +000.02 000.0 81 0.999 +000.11 
Be 1.002 -000.02 000.0 82 0.999 +000. 12 
33 1.002 -000.06 359.9 83 1.000 +000. 13 
34 1.001 -000.09 359.9 84 1.001 +000.11 
35 1.001 2000 sie 359.9 85 1.002 +000.09 
36 1.000 000.43 359.9 86 1.002 +000. 06 
33 0.999 00012 359.9 87 1.002 +000.02 
38 0.999 -000.10 359.9 88 1,002 -000.03 
39 0.998 -000.07 359.9 89 1.002 ~000.07 
1940 0.998 -000.03 000.0 1990 1.001 -000.10 
41 0.998 +000.01 000.0 91 1.001 =000;.92 
42 0.998 +000.05 000.0 92 1.000 ~000.13 
43 0.998 +000.09 000.1 93 0.999 -000.12 
44 0.999 +000.11 000.1 94 0.999 -000.10 
45 1.000 +000.12 000.1 95 0.998 -000.07 
46 1.000 +000.12 000.1 96 0.998 -000.03 
47 1.001 +000.11 000.1 97 0.998 +000.01 
48 1.002 +000.08 000.1 98 0.998 +000. 06 
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11. FIGURE CAPTIONS 


Frontispiece. Place names in Canada on a Miller projection. 


Fig. I. Relation between the local constituent A cos (ot-a) and the 
corresponding component of the tidal force A' cos V(t). 


Fig. 2. Profile of the surface of a fluid in a canal and corresponding 
horizontal velocities at various times when it is disturbed by (a) a 
travelling wave and (b) a standing wave. 


Fig. 3. Cotidal charts in a canal half a wavelength long (one-dimensional 
motion). Continuous line - cophase contour. Dashed line - coamplitude 
contour. 


a) Travelling wave: vertical displacement Z. 

b) Standing wave, no friction: vertical displacement Z and 
horizontal velocity u. 

c) Standing wave affected by friction: Z and u. 


Fig. 4. Schematic cotidal charts for a standing Kelvin wave (no friction); 
they are the rotary counterpart of chart (b) in Fig. 3. The nodal line has 
been transformed into a point (amphidromy) by the rotation and the cophase 
lines now radiate from it. The node of the current (at X]) corresponds to 
the antinode of the vertical displacement and vice versa. 


Fig. 5. Cotidal chart of the vertical displacement Z for a standing Kelvin 
wave corresponding to the frequency of M2 in a basin at latitude 45°N, width 
200 km, depth 100 m, closed at x = 0 and of indefinite_length. (a) Without 
friction. (b) With weak linearized friction, r = 10-°/sec. (c) With 
moderate friction, r = 5 x 10-9/sec. (d) With strong friction, r = 107'*/sec. 


Fig. 6. Set of predicted levels for April 1978 at various stations to 
illustrate some actual tides corresponding to the indicated patterns. The 
stations chosen cover patterns from .4SD (almost pure semi diurnal) to 2.0D 
(almost pure diurnal). (a) 391 SD. Woods Harbour. Good semidiurnal tide 
with a slight diurnal inequality. (b) .483 SD. Halifax. Quite Similar to 
Woods Harbour with a more marked diurnal inequality and more intense change 
from springs to neaps. (c) .608 SD. Riviére du Loup. Still a good semi- 
diurnal tide with an even more marked change from springs to neaps. 

(d) 1.004 N. Holy Rood. We enter the normal category where the es 
of the constituents is very similar to that of the tidal potential. (e eee 
Rivi@re au Renard. We enter the M (mixed) category where Be, SAT Pit O° 
of the spectrum of the constituents becomes larger than that in the tida ‘= 
potential. The diurnal inequality is quite marked during the two ae: 
maxima in the lunar declination which occur around hour 288 and hour ne 

the month. (f) .306 M. Pictou. (g) .501 M. Seattle, Washington, 


ae 


(h) .602 M. Lund. (i) .802 M. Gabriola Pass. (j) 1.003 D. Cowichan. 
We enter the D (diurnal) category where the probability increases that there 
will occur a single tide over a span of 24 hours. (k) 1.205 D. Saanichton 
Bay.(2) 1.401 D.. EsquimaltsAi(mjgnim603sb2 .GaissteyPoints Ihbsrstation 
lies to the southwest of the semidiurnal point of amphidromy which exists in 
Northumberland Strait. (n) 2.021 D. Gorge Victoria. Apparently the record 
with the most diurnal character available in our files. Both the diurnal 

and semidiurnal bands of the mixed tide entering from Victoria are damped 
while it propagates into the Gorge but the semidiurnal band is damped rapidly. 


Fig. 7. Observed amplitude ratio and difference in the Greenwich phase lag 
between the pairs K2 and So, P] and Kj. It is presented as a vector of 
length r (amplitude ratio) and inclined at an angle o (difference of phase 
lag). These values are obtained at stations where the constituents are 
separable; they are used afterwards for inferring K2 from So and P} from Kj] 
at stations where they cannot be separated. Since the observed values depend 
on the geographical location, the measured vectors are presented in groups 
for the west coast, east coast and arctic. The tendency of the vector to 
bunch around a common value indicates that a mean value exists for the vector 
and the vector mean is presented in each case. The actually observed r, 9 
values are shown as rays in the larger diagram. Their vector average is 
shown in the smaller diagram as a single ray whose amplitude (r) and orien- 
tation (¢, in degrees) are written out explicitly as "mean amplitude" and 
"mean angle". 


Fig. 8. Cotidal charts for the constituents M2, S?, N2, K] and 0} in waters 
adjacent to Canada. A Miller projection is used for the geographic back- 
ground. The solid lines are contours of constant Greenwich phase lag; the 
dashed lines are contours of constant amplitude and the units degrees and 
centimeters. The contours are marked as 120° or 30 to indicate a Greenwich 
phase lag of 120° or a contour of 30 cm. The time zone used for all the 
maps is Greenwich (Z = 0). The spacing of the contours is not necessarily 
regular. Question marks indicate areas where it was not possible to con- 
jecture the type of motion present either because of a shortage of data or 
of their poor quality. The signs + or - help visualize the direction of 
increasing or decreasing values in areas where this is not evident. (a) M2 
Cy Byam O08 i Che 0 le Pah 9) 8 


Fig. 9. Tide patterns in waters adjacent to Canada. 


Fig. 10. Cotidal charts for the constituents M2, S92, No, K] and 07 in waters 
adjacent to the northern and southern portions of British Columbia. (a) Moe 
(b) sae (ChoaNo AdeckT CejesdT: 


Fig. 11. Patterns of the tide in waters adjacent to the northern and 
southern portions of British Columbia. (a) Northern portion. (b) Southern 
portion. 


Fig. 12. Cotidal charts for the constituent M2, S2, No, Kj and 0] around 
Newfoundland and off Cape Breton to the estuary of the St. Lawrence River. 
(a) Ho (bh 52 (er Na Ad At ee ue 


Fig. 13. Cotidal charts for the constituent Mo, S2, N2, K] and 0] in the 
upper St. Lawrence River. (a) M2 (b) Sa (c) No (d) Ky (e) 0}. 
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Fig. 14. Creation of a semimonthly tide in the upper reaches of the St. 
Lawrence River caused by the succession of spring and neap tides downstream 
(from Godin 1979). 


Fig. 15. Tidal pattern around the Gulf of St. Lawrence and the upper 
estuary of the St. Lawrence River. (a) Gulf of St. Lawrence (b) Upper 
estuary. 


Fig. 16. Cotidal charts for the constituents Mo, So, No, Ky and 0] in the 
Bay of Fundy. (a) Mo (b) Sp (c) No (d) Ky a 0}. 


